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Chimera states consisting of domains of coherently and incoherently oscillating nonlocally-coupled 
phase oscillators in systems with spatial inhomogeneity are studied. The inhomogeneity is intro¬ 
duced through the dependence of the oscillator frequency on its location. Two types of spatial 
inhomogeneity, localized and spatially periodic, are considered and their effects on the existence 
and properties of multi-cluster and traveling chimera states are explored. The inhomogeneity is 
found to break up splay states, to pin the chimera states to specific locations and to trap traveling 
chimeras. Many of these states can be studied by constructing an evolution equation for a complex 
order parameter. Solutions of this equation are in good agreement with the results of numerical 
simulations. 

PACS numbers: May be entered using the \pacs{#l} command. 


I. INTRODUCTION 

Networks of coupled oscillators have been extensively 
studied for many years, owing to their wide applica¬ 
bility in physics, chemistry, and biology. As examples 
we mention laser arrw^ Josephson junctions, popula¬ 
tions of fireflies, etc. [ll-Q- The phase-only models have 
proved to provide useful models for systems with weak 
coupling. The best known model of this type is the Ku- 
ramoto model in which the oscillators are described by 
phase variables ^ and coupled to others through a sinu¬ 
soidal function [3-Q. These models exhibit a transition 
to collective synchronization as the coupling strength in¬ 
creases, a process that has been described as a phase 
transition. A general form of these systems is as follows: 

(to k ^ 

-^ = Wi - — ^ Gy sin(0, - Oj + a). (1) 

j=t 

Here tOi is the natural frequency of oscillator i, Gij rep¬ 
resents the coupling between oscillators i and j, a is a 
phase lag and k is the overall coupling strength. A gen¬ 
eral treatment of this system is not easy, and two types 
of simplifications are commonly used. One of them is to 
assume global coupling among the oscillators, and assign 
the natural frequencies uii randomly and independently 
from some prespecified distribution (MU. The second 
tractable case arises when all the oscillators are assumed 
to be identical, and the coupling G, y is taken to be lo¬ 
cal, eg., nearest-neighbor coupling [Q. The intermediate 
case of nonlocal coupling is harder but the phenomena 
described by the resulting model are much richer. 

In this paper we suppose that the oscillators are ar¬ 
ranged on a ring. In the continuum limit A —>■ oo the 
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system is described by the nonlocal equation 

00 f 

— = Uj(x) - I G{x - y) sin[6'(a:, t) - 9{y, t) + a] dy (2) 

for the phase distribution 9{x,t). When the oscillators 
are identical (w is a constant) and G{x) = ^ exp(— k|x|) 
this system admits a new type of state in which a frac¬ 
tion of the oscillators oscillate coherently (i.e., in phase) 
while the phases of the remaining oscillators remain inco¬ 
herent [l3j| . In this paper we think of this state, nowadays 
called a chimera state 0 , as a localized structure em¬ 
bedded in a “turbulent” background. Subsequent studies 
of this unexpected state with different coupling functions 
G{x) have identified a variety of different one-cluster and 
multi-cluster chimera states [l5l - l2^ , consisting of clusters 
or groups of adjacent oscillators oscillating in phase with 
a common frequency. The clusters are almost stationary 
in space, although their position (and width) fluctuates 
under the influence of the incoherent oscillators on either 
side. Recently, a new type of chimera state has been 
discovered, a traveling chimera state [1^. In this state 
the leading edge plays the role of a synchronization front, 
which kicks oscillators into synchrony with the oscillators 
behind it, while the trailing front kicks oscillators out of 
synchrony; these two fronts travel with the same speed, 
forming a bound state. 

It is natural to ask whether these states persist in 
the presence of spatial inhomogeneity, in particular, in 
the presence of spatial inhomogeneity in the natural fre¬ 
quency distribution (w = w(a:)). In [T^ Laing demon¬ 
strated the robustness of the chimera state with re¬ 
spect to inhomogeneity in a two population model, which 
can be considered to be the simplest model exhibiting 
chimera states [T^. However, in this model there is no 
spatial structure to either population. Consequently, we 
focus in this paper on a ring of adjacent oscillators with a 
prescribed but nonuniform frequency distribution ui(x), 
where the continuous variable x (—tt < x < tt) repre¬ 
sents position along the ring. To be specific, two classes 
of inhomogeneity are considered, a bump inhomogeneity 
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uj(x) = Wo exp(—K|a;|), k > 0, and a periodic inhomo¬ 
geneity ijj{x) = Wn cosjlx), where I is a positive integer. 
In each case we follow and study the coupling func¬ 
tions 

= cos{nx), 

G^^x) = cos(nx) + cos[(n -I- l)a;], 

where n is an arbitrary positive integer. These choices 
are motivated by biological systems in which coupling be¬ 
tween nearby oscillators is often attractive while that be¬ 
tween distant oscillators may be repelling . There are 
several advantages to the use of these two types of cou¬ 
pling. The first is that these couplings allow us to obtain 
chimera states with random initial conditions. The sec¬ 
ond is that we can identify a large variety of new states for 
suitable parameter values, including (a) splay states, (b) 
stationary multi-cluster states with evenly distributed co¬ 
herent clusters, (c) stationary multi-cluster states with 
unevenly distributed clusters, (d) a fully coherent state 
traveling with a constant speed, and (e) a single cluster 
traveling chimera state [ 2 ^ . 

In this paper, we analyze the effect of the two types 
of inhomogeneity in w on each of these states and de¬ 
scribe the results in terms of the parameters ujq and k or 
I describing the strength and length scale of the inhomo¬ 
geneity, while varying the parameter l3 = ^—a represent¬ 
ing the phase lag a. In Section |TT1 we briefly review the 
notion of a local order parameter for studying chimera 
states and introduce the self-consistency equation for this 
quantity. In Sections Ell El and El we investigate, re¬ 
spectively, the effect of inhomogeneity on rotating states 
(including splay states and stationary chimera states), 
traveling coherent states and the traveling chimera state. 
We conclude in Section IVII with a brief summary of the 
results and directions for future research. 


II. EFFECTIVE EQUATION 

Equation ([2|) is widely used in studies of chimera states. 
An equivalent description can be obtained by construct¬ 
ing an equation for the local order parameter z{x,t) de¬ 
fined as the local spatial average of exp[f0(a;, t)], 

1 

z{x,t)= lim - / (3) 

( 5 ^ 0 + 0 J_s/2 

The evolution equation for z then takes the form [28l - [^ 

zt = iuj{x)z + i - e^z^Z*) , (4) 

where Z{x,t) = A'[z](x,t) and if is a compact linear 
operator defined via the relation 

K[u]{x,t)= f G{x -y)u{y,t)dy. (5) 

J —TT 

A derivation of Eq. @ based on the Ott-Antonsen 
Ansatz jEl is given in the Appendix. Equation (jT]) can 


also be obtained directly from Eq. m using the change 
of variable 


z{x,t) = exp[i6{x,t)]. (6) 


An important class of solutions of Eq. © consists of 
stationary rotating solutions, i.e., states of the form 

z(a:, t) = z(a;)e“*^*, (7) 


whose common frequency O satisfies the nonlinear eigen¬ 
value relation 


i\p. + oj{x)\z + - e ^°'Z{x) — z^e‘^Z*{x) =0. (8) 


Here z{x) describes the spatial profile of the rotating so¬ 
lution and Z = K[z]. 

Solving Eq. ([8]) as a quadratic equation in z we obtain 


Z*(x) 


e^^Z(x) 

Q + u(x) + y.{x) 


(9) 


The function y is chosen to be [(Q -I- — \Z\'^YG when 

|Q-|-w| > \Z\ and i[\Z\^ — when |Q-|-a;| < \Z\. 

This choice is dictated by stability considerations, and in 
particular the requirement that the essential spectrum 
of the linearization about the rotating solution is either 
stable or neutrally stable [ 2 ^ . The coherent (incoherent) 
region corresponds to the subdomain of (—tt, tt] where 
|Q -I- oj{x)\ falls below (above) \Z{x)\. Substitution of 
expression (|9]) into the definition of Z{x) now leads to 
the self-consistency relation 

= + ( 10 ) 

Here the bracket (•) is defined as the integral over the 
interval [—7r,7r], i.e., 

{u)= f u{y)dy. (11) 


In the following we write 

Z{x) = R{x)e^^^^\ (12) 


and refer to R{x) and ©(x) as the amplitude and phase 
of the complex order parameter Z{x). 


III. STATIONARY ROTATING SOLUTIONS 

As shown in [ 2 ^, Eq. (I2|) with constant natural fre¬ 
quency uj exhibits both stationary rotating solutions 
(splay states and stationary chimera states) and traveling 
solutions (traveling coherent states and traveling chimera 
states) for suitable coupling functions G{x). In this sec¬ 
tion, we investigate the effect of spatial inhomogeneity 







3 


in w (i.e., uj = w(x)) on the splay states and on station¬ 
ary chimera states, focusing on the case G^\x) = cosnx 
studied in [2^. We find that when the inhomogeneity 
is sufficiently weak, the above solutions persist. How¬ 
ever, as the magnitude of the inhomogeneity increases, 
new types of solutions are born. The origin and spatial 
structure of these new states can be understood with the 
help of the self-consistency relation (US. Since Eq. (nni) 
is invariant under the transformation 2; ze''’^ with ^ an 

arbitrary real constant, the local order parameter Z(x) 
for the coupling function G{x) = cosnx can be written 
as 


Z{x) = a cos nx + b sin nx, 


(13) 



where a is positive and b = br + ibi with br and bi both 
real. Substituting the Ansatz (fT51) into Eq. (flUl) we obtain 
the following pair of integral-algebraic equations 

cos{ny){n{y) - ^i{y)) 
a cos ny + b* sin ny 

sin{ny)(Sl{y) - y,(y)) 
a cos ny + b* sin ny 

where y{y) = — |acos(n2/) -|- bsm{ny)\'^y/^ and 

r2(?/) = + w(?/). These equations may also be written 

in the more convenient form 

(^(y)-M(y)) =e"*'^(a^ + |6p), (16) 

sin{2ny) — 2ab cos{2ny) 

^(y)+A»(y) 

In the following we consider two choices for the inho¬ 
mogeneity w(x), a bump w(x) = wq exp (—«;|x|), k > 0, 
and a periodic inhomogeneity uj{x) = ujocos{lx) where 
I is a positive integer. The resulting equations possess 
an important and useful symmetry, x —x, b —b. 
Moreover, for uj(x) satisfying a;(x -I- 2t:/ 1) = uj(x) with I 
an integer, a solution Z(x) of the self-consistency relation 
implies that Z{x + 2nnv/l) is also a solution. Here m is 
an arbitrary integer. 

In the following we use Eqs. (USD and (HZl repeatedly 
to study the changes in both the splay states and the 
stationary chimera states as the magnitude wq > 0 of the 
inhomogeneity increases, and compare the resulting pre¬ 
dictions with numerical simulation of = 512 oscillators 
evenly distributed in [—7r,7r]. 






A. Bump inhomogeneity: a; = cao exp ( —k|x|) 

We now consider the case where uj(x) has a bump de¬ 
fect at X = 0. For the spatial profile of the defect we 
pick a;(x) = ujq exp (— k|x|), where wq > 0 and k > 0 are 
parameters that can be varied. 


FIG. 1. The phase distribution 9{x) for splay states observed 
with G{x) = cos(x), oj{x) = wo exp(—2|x|) and b = 0.05. (a) 
Wo = 0.006. (b) wo = 0.02. (c) wo = 0.05. (d) wo = 0.1. The 
states travel to the right (H > 0). 


1. Effect on splay states 

When G{x) = cos(nx) and w is a constant, Eq. (I5|) 
exhibits so-called splay state solutions with 6{x, t) = qx — 
fit, where fl is the overall rotation frequency and q is 
an integer called the twist number. The phase in this 
type of solution drifts with speed c = fl/q to the right 
but the order parameter is stationary. Consequently we 
think of the splay states as a stationary rotating states. 
Linear stability analysis for G(x) = cosnx shows that 
the splay state is stable when |g| = n [1^, a result that is 
easily confirmed in simulations starting from randomly 
distributed initial phases. In the following we consider 
the case n = 1 and simulate N = 512 oscillators evenly 
distributed in [—7r,7r] for different values of wq and k. 

Figure [T] shows how the splay solutions change as ojq 
increases. When ojq becomes nonzero but remains small, 
the splay states persist but their phase 9{x,t) no longer 
varies uniformly in space (Fig. [T](a)). Instead 

6{x, t) = 4>{x) — fit, (18) 

where ^(x) is a continuous function of x with (piir) — 
= 2'Kq. We refer to this type of state as a near¬ 
splay state. As ojq becomes larger an incoherent region 
appears in the vicinity of x = 0, with width that in¬ 
creases with increasing wq (Figs. [IKb,c)). We refer to 
this type of state as a chimera splay state. As ojq in¬ 
creases further and exceeds a second threshold, a new 
region of incoherence is born (Fig. [T](d)). Figure [2] pro¬ 
vides additional information about the partially coherent 
near-splay state in Fig. [Hb). The figure shows the real 
order parameters i?(x) and 0(x) together with 9t, the lo¬ 
cal rotation frequency averaged over a long time interval 
(Fig. HKd)), and reveals that in the coherent region the 
oscillation frequencies are identical (with 9t = —H), with 
an abrupt but continuous change in the frequency dis- 
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FIG. 2. (a) A snapshot of the phase distribution 9{x,t) for 
G{x) = cos{x) and uj{x) = 0.02exp( —2|a:|). (b) The corre¬ 
sponding R{x). (c) The corresponding 0(x). (d) The oscilla¬ 
tor frequency dt averaged over the time interval 0 < t < 1000 
(solid line). In the coherent region coincides with the global 
oscillation frequency —fl (open circles). All simulations are 
done with /3 = 0.05 and N = 512. 


tribution within the incoherent region. The figures also 
reveals that as wq increases the incoherent region devel¬ 
ops a stronger and stronger asymmetry with respect to 
X = 0, the bump maximum. This is a consequence of the 
asymmetry introduced by the direction of travel, i.e., the 
sign of the frequency fl in Eq. m as discussed further 
below. 

These states and the transitions between them can 
be explained within the framework of the self-consistent 
analysis. To compute the solution branches and the tran¬ 
sition thresholds, we numerically continue solutions of 
Eqs. dm) and (El) with respect to the parameter ujq. 
When ujQ = 0, the splay state (with positive slope) cor¬ 
responds to a = TT, br = 0, bi = TT, with fl = tt. The 
order parameter is therefore i?exp(i0) = 7rexp(±*a;). 
We use this splay state as the starting point for continu¬ 
ation. As Wo increases, a region of incoherence develops 
in the phase pattern in the vicinity of a: = 0. From the 
point of view of the self-consistent analysis, the incoher¬ 
ent region corresponds to the region where the natural 
frequency exceeds the amplitude R{x) of the complex or¬ 
der parameter. The boundaries between the coherent and 
incoherent oscillators are thus determined by the relation 
il + u}{x) = |acos(x) -I- {br + ibi) sin(a;)|. For wq = 0.02, 
the left and right boundaries are thus xi ~ —0.4026 and 
Xr ~ 0.2724, respectively. These predictions are in good 
agreement with the values measured in direct numeri¬ 
cal simulation (Fig. [2](a)). Figure [31(a) shows the overall 
frequency fl obtained by numerical continuation of the 
solution of Eqs. and ED in the parameter wq , while 
Figs. |3Db,c) show the corresponding results for the frac¬ 
tion e of the domain occupied by the coherent oscillators 
and the extent of the (first) region of incoherence, i.e., the 
interval xi < x < Xr, also as functions of wq. From the 


(a) 



(b) 



(c) 



FIG. 3. (a) The overall frequency SI, (b) the fraction e of 
the domain occupied by the coherent oscillators, and (c) the 
width of the incoherent region xi < x < Xr, all as functions 
of the parameter wq. The calculation is for k = 2, /3 = 0.05 
and N = 512. 


figure we can see a clear transition at wq Ri 0.0063 from 
a single domain-filling coherent state to a “splay state 
with one incoherent cluster” in xi < x < Xr, followed 
by a subsequent transition at ojq w 0.065 from this state 
to a “splay state with two incoherent clusters”. These 
transitions occur when the profiles of $7 -|-w(a:) and R(x) 
touch as Wq increases and these points of tangency there¬ 
fore correspond to the locations where coherence is first 
lost. Figure ID shows that tangencies between fl -|- w(a:) 
and R(x) occur when wq ss 0.0063 and 0.065, imply¬ 
ing that intervals of incoherent oscillators appear first at 
a; = 0 (i.e., the bump maximum) and subsequently at 
X Ri 1.83, as Wo increases. These predictions are in ex¬ 
cellent agreement with the direct numerical simulations 
shown in Fig.[T] Moreover, the critical values of wq pre¬ 
dicted by the self-consistency analysis are fully consistent 
with the simulation results when /3 is increased quasi- 
statically (not shown). In each case we repeated the 
simulations for decreasing /3 but found no evidence of 
hysteresis in these transitions. 
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FIG. 4. Comparison of n+ti;(a;) and R{x) at the critical values 
tjo for the appearance of new regions of incoherence around 
(a) a; = 0 for ljq ~ 0.0063 and (b) x = 1.83 for luq = 0.065. 
The calculation is for ac = 2, /? = 0.05 and N = 512. 

2. Effect on stationary chimera states 

Chimera states with 2n evenly distributed coherent 
clusters are readily observed when G{x) = cos{nx) and 
w is a constant. These states persist when a bump is in¬ 
troduced into the frequency distribution uj{x). Figure O 
shows the phase distribution and the corresponding local 
order parameters R{x) and 0(a;) for G{x) = cos{x) when 
uj{x) = 0.1 exp(—2|a;|). The figure shows that the two 
clusters persist, but are now always located near x = —^ 
and j. This is a consequence of the fact that the presence 
of the bump breaks the translation invariance of the sys¬ 
tem. Figures ElJbjC) show that the local order parameter 
Z(x) has the symmetry Z{—x) = —Z(x). This symme¬ 
try implies Z should take the form i?exp (10) = &sin(a::), 
where b = br is real. The corresponding self-consistency 
equation takes the form 

-I- uj{y) — Y^(r2 -I- uj{y)Y — b"^ sin^ . 

(19) 

The result of numerical continuation of the solutions of 
Eq. (fT^ are shown in Fig.[6l The 2-cluster chimera state 
persists to large values of ujq, with the size of the coher¬ 
ent clusters largely insensitive to the value of ujq. This 
prediction has been corroborated using direct simulation 
of Eq. ([T|) with N = 512 oscillators. 

When uj is constant, finite size effects cause the phase 
pattern to fluctuate in location. In [^ . we demonstrate 
that this fluctuation is well modeled by Brownian motion 
in which the variance is proportional to t, even though the 
original system is strictly deterministic. As mentioned 


(c) 


2 

® 0 ^ 

-2 

-2 0 2 
X 

FIG. 5. (a) A snapshot of the phase distribution 6(x,t) in 
a 2-cluster chimera state for G{x) = cos(x) and uj(x) = 
0.1 exp(—2|x|). (b) The corresponding order parameter R(x). 
(c) The corresponding order parameter Q>{x). Note that the 
oscillators in the two clusters oscillate with the same fre¬ 
quency but TT out of phase. The calculation is done with 
K = 2, /3 = 0.05 and N = 512. 


above, when ijj{x) is spatially dependent, the translation 
symmetry is broken and the coherent cluster has a pre¬ 
ferred location. Figure [5] suggests that local maxima of 
the order parameter R can be used to specify the location 
of coherent clusters. Consequently we plot in Fig. [7] the 
position XQ{t) of the right coherent cluster as a function 
of time for three different values of the parameter k. We 
see that the inhomogeneity pins the coherent cluster to a 
particular location, and that the cluster position executes 
apparently random oscillations about this preferred lo¬ 
cation, whose amplitude increases with increasing k, i.e., 
with decreasing width of the bump. Figure |S] shows the 
standard deviation of the position XQ(t) of the coherent 
cluster as a function of the parameter k. 

The behavior shown in Figs. [7] and [5] can be mod¬ 
eled using an Ornstein-Uhlenbeck process, i.e., a linear 
stochastic ordinary differential equation of the form 

dxo = X{fi — xo)dt + adWt, (20) 
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FIG. 6. The dependence of (a) Q and the order parameter 
amplitude b on loq. (b) The fraction e of coherent oscillators 
as a function of uq. The calculation is done with k = 2, 
j3 = 0.05 and N = 512. 


;o 

CO 


0 . 2 |-^-- 

I ° 

0.1 

o 

I ° 

0 

I o 

°0 ^ 5 10 


K 


FIG. 8. The standard deviation of the position xci{t) as a 
function of the parameter k for cjq = 0.1, /3 = 0.05 and N = 


512. 
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FIG. 7. The position xo{i) of the coherent cluster as a func¬ 
tion of time when cjq = 0.1, /3 = 0.05, N = 512 and (a) 
K = 10, (b) K = 6 and (c) k = 2. 


FIG. 9. The histogram of the residues e = xo{t + At) — 
Axo{t) — B is well approximated by a normal distribution 
when At = 100. 

where A represents the strength of the attraction to the 
preferred location /r, and cr indicates the strength of the 
noise. Models of this type are expected to apply on an ap¬ 
propriate timescale only: the time increment At between 
successive steps of the stochastic process must be large 
enough that the position of the cluster can be thought of 
as the result of a large number of pseudo-random events 
and hence normally distributed, but not so large that 
nonlinear effects become significant. Figure |9] reveals 
that for an appropriate interval of At the fluctuations 
xo{t + At) — Axo{t) — B are indeed normally distributed, 
thereby providing support for the applicability of Eq. (1^ 
to the present system. 

Equation ((20l) has the solution 

xo(t -f At) = Axoit) B CNq i, (21) 

where A = B = p,(l — A) and C = 

To fit the parameters to the data in Eig. [3 we notice 
that the relationship between consecutive observations 
xq is linear with an i.i.d. error term CAo.i, where Aq.i 
denotes the normal distribution with zero mean and unit 
variance (Eig. [5]) and C is a constant. A least-squares fit 
to the data {xo{t),xo{t + At)) gives the parameters A, /r 
and a. We And that for ujq = 0.1, k = 2 and P = 0.05, 
the choice At = 100 works well and yields the empirical 
model parameters A ~ 0.016, fv —1.57 and a ~ 0.02, 
a result that is in good agreement with the simulation 
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(a) 






FIG. 10. A near-splay state for (a) oj{x) = 0.1cos(2x). 
Chimera splay states for (b) lo{x) — 0.2cos(2a;), (c) lo{x) = 
0.3cos(2a;) and (d) uj{x) = 0.4cos(2a;). In all cases /3 = 0.05 
and N = 512. 
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results for At = 2 summarized in Fig. [5] 


B. Periodic io{x) 

In this section we consider the case uj{x) = ujocos{lx). 
When Wo is small, the states present for wq = 0 persist, 
but with increasing wq one finds a variety of intricate 
dynamical behavior. In the following we set G{x) = cos a; 
and j3 = j — a = 0.05, with ujq and I as parameters to 
be varied. 



FIG. 11. Chimera splay states for G{x) = cos(a;) and (a) 
oj{x) = 0.2cos(3a;) (3-cluster state), (b) uj{x) = 0.2cos(4a;) 
(4-cluster state), and (c) w(x) = 0.2cos(5a:) (5-cluster state). 
In all cases (3 = 0.05 and N = 512. 


1. Splay states, near-splay states and chimera splay states 

When oj{x) is a constant, splay states are observed in 
which the phase 9{x,t) varies linearly with x. When wq 
is nonzero but small, the splay states persist as the near¬ 
splay states described by Eq. (HID; Fig. fTOT a) shows an 
example of such a state when u}{x) = 0.1cos(2x). As 
Wo becomes larger, incoherent regions appear and these 
increase in width as wo increases further (e.g., Fig. llOf b)- 
(d)). We refer to this type of state as a chimera splay 
state, as in the bump inhomogeneity case. Figures [TUI 
show the phase distribution obtained for different values 
of wq. For these solutions, the slope of the coherent region 
is no longer constant but the oscillators rotate with a 
constant overall frequency VL. This type of solution is also 
observed for other values of 1. Figure [IT] shows examples 
of chimera splay states for Z = 3, 4 and 5, with I coherent 
clusters in each case. 

To compute the solution branches and identify thresh¬ 
olds for additional transitions, we continue the solutions 
of Eqs. (fTUl) and (fTTI) with respect to the parameter wq. 
Figures (TUI and (TUI show the dependence of and of the 
coherent fraction e on wq when I = 1 and 2, respectively. 
The figures indicate that the coherent fraction e falls be¬ 
low 1 at Wo ss 0.0075 {I = 1), wq ~ 0.16 {I = 2), and 


Wq w 0.0048 (Z = 3, 4 and 5). These values coincide 
with the parameter values at which an incoherent region 
emerges. In addition. Fig. ITUTb) reveals the presence of 
a second transition, at wq sa 0.13, corresponding to the 
emergence of a 2-cluster state from a 1-cluster state. Fig¬ 
ures M and M show the profiles of -|- w(x) and R{x) 
at the critical values of ojq at which incoherent regions 
emerge. 

To understand the Z-cluster chimera states reported in 
Fig. HU we computed the local order parameter Z(x) 
and found that i?(x) is approximately constant while 
the phase 0(x) varies at a constant rate, suggesting the 
Ansatz Z(x) = ae“. With this Ansatz, Eqs. (fTU|) and 
(HZD reduce to a single equation, 

27rQ - + Wo cosily))^ - (22) 

for all positive integers I > 3. Figure (TUI shows the result 
of numerical continuation of a solution for Z = 3. The 
figure reveals no further transitions, indicating that the 
I = 3 chimera state persists to large values of Wq. Nu¬ 
merical simulation shows that these states (Z = 3, 4, 5) 
are stable and persist up to wq = 1. 












FIG. 12. The dependence of (a) fl and (b) the fraction e of 
the domain occupied by the coherent oscillators on ujq, with 
uj{x) = ojo cos(x) and P = 0.05. 


(a) 



(b) 



FIG. 13. The dependence of (a) and (b) the fraction e of 
the domain occupied by the coherent oscillators on ujq, with 
(jj{x) = ioo cos(2a:) and /3 = 0.05. 


FIG. 14. Gomparison of + oj{x) and R{x) at the critical 
values ujQ for the appearance of a new region of incoherence 
around (a) a; = 0 for ojq = 0.0075 and (b) x = 2.8 for uo = 
0.13. Parameters: I = 1 and P = 0.05. 



X 


FIG. 15. Gomparison of + uj{x) and R{x) at the critical 
value oiQ ~ 0.16. Parameters: I = 2 and P = 0.05. 


1-cluster state we use Fig.lYTTbl to conclude that Z{x) is 
of the form acos(a:), and hence that Eqs. (fTCl) and (fT71) 
reduce to the single equation 

27rf2 — + oj{y)Y — cos^ (23) 


2. 1-cluster chimera states 

Figure[T7]shows an example of a 1-cluster chimera state 
when u}{x) = 0.1cos(a::); 1-cluster chimera states of this 
type have not thus far been reported for G{x) = cos(a;), 
Wq = 0, where computations always result in 2-cluster 
states [^. To understand the origin of this unexpected 


Figure [18] shows the result of numerical continuation of 
the solution of this equation as a function of ojq- The 
figure reveals a transition at ujq = 0.0463 iFigs. fTST a.b')'). 
When ujo < 0.0463, the solution is a 2-cluster chimera; 
as ujq increases through ujq ~ 0.0463 the coherent region 
around x = 0 disappears, leaving a single cluster chimera 
state (Fig. [THKc)). 
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(a) 



(b) 



FIG. 16. Dependence of (a) Q and (b) the coherent fraction 
e on ajo, with uj{x) = loq cos(3a;) and P = 0.05. 


(a) 



X 


FIG. 17. (a) A snapshot of the phase distribution 6(x, t) 

in a 1-cluster chimera state for G(x) = cos( 2 ;) and lij(x) = 
0.1 cos (a;), (b) The corresponding order parameter R{x). (c) 
The corresponding order parameter 0(a;). Figure (a) shows 
that the presence of a nearly coherent region near x = 0 with 
oscillators that oscillate tt out of phase with the coherent clus¬ 
ter. The calculation is done with /3 = 0.05 and N = 512. 


(a) 



(c) 



FIG. 18. (a) The quantities a and Q, and (b) the coherent 
fraction e, all as functions of wo for the chimera state shown 
in Fig. 1171 (c) D -1- oj{x) and R{x) as functions of x when 
LUO = 0.0463. 


3. 2-cluster chimera states 

For G{x) = cos{x) and constant w simulations always 
evolve into either a 2-cluster chimera state or a splay state 
[2^. It is expected that when wq is small, the 2-cluster 
chimera is not destroyed. Figures [TST a.bl gives examples 
of this type of state with I = 2 and I = 3 inhomogeneities, 
respectively. In both cases the clusters are located at 
specific locations selected by the inhomogeneity. 

The order parameter profile with I = 3 correspond¬ 
ing to the numerical solution in Fig. ITW bl has the same 
symmetry properties as that in Figs. [Sljb,c). It follows 
that we may set a = 0,6 = 6r>0in the order pa¬ 
rameter representation 0, resulting once again in the 
self-consistency relation 0. We have continued the so¬ 
lutions of this relation for I = 3 using the simulation in 
Fig. Iiyf bl with Wo = 0.1 to initialize continuation in wq. 
Figure [20] shows the result of numerical continuation of 
the order parameter for this state. No further transitions 
are revealed. There are three preferred locations for the 
coherent clusters, related by the translation symmetry 
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FIG. 19. 2-cluster chimera states for G(x) = cos(a;) and (a) 
iij(x) — 0.1cos(2a;) and (b) u!(x) = 0.1cos(3a:). In both cases 
/3 = 0.05 and N = 512. 


(a) 



(b) 



“o 


FIG. 20. The dependence of (a) SI (solid line) and b (dashed 
line), and (b) the coherent fraction e on ujo when I = 3 and 
/3 = 0.05. 
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FIG. 21. The three possible locations of the 2-cluster chimera 
state when uj{x) = 0.1cos(3a:). The simulation is done with 
P = 0.05 and N = 512. 

(a) (b) 




FIG. 22. (a) A snapshot of the phase pattern in a traveling 
coherent state when oj is constant, (b) The position xq of 
this state as a function of time. The simulation is done for 
G{x) = cos(x) -I- cos(2a;) with P = 0.75 and N = 512. 


X ^ x + 2^, m = 0,1, 2, as shown in Fig. HT] 

IV. TRAVELING COHERENT SOLUTIONS 

We now turn to states with a spatially structured or¬ 
der parameter undergoing translation. With G(x) = 
cos(a::) + cos(2a::) and constant uj the system ([5]) exhibits 
a fully coherent but non-splay state that travels with 
a constant speed c(/3) when 0.646 3 3 3 0.7644 [^ . 
Figure E^ a) shows a snapshot of this traveling coherent 
state while Fig. E^ b) shows its position xq as a func¬ 
tion of time. The speed of travel is constant and can be 


obtained by solving a nonlinear eigenvalue problem [^ . 
In this section, we investigate how the inhomogeneities 
uj{x) = Wo exp (—K|g:|) and u}{x) = ujocos(lx) affect the 
dynamical behavior of this state. 

A. Bump inhomogeneity: w(x) = wo exp (—k|x|) 

In this section we perform two types of numerical ex¬ 
periments. In the first, we fix /3 and start with the trav¬ 
eling coherent solution for wq = 0. We then gradually 
increase wq gradually at fixed values of k and p. In the 
second we fix wq > 0, k > 0, and vary p. 
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For each 0.646 < /3 < 0.7644 we find that the coher¬ 
ent state continues to travel, albeit nonuniformly, until 
Wo reaches a threshold value that depends of the values 
of /3 and k. The case (3 = 0.75 provides an example. 
Figure ISST al shows the traveling coherent state in the 
homogeneous case (wo = 0) while Fig. [23fbl shows the 
corresponding state in the presence of a frequency bump 
Cij(x) = Wo exp (—K|a:|) with wq = 0.04 and k = 2. In 
this case the presence of inhomogeneity leads first to a 
periodic fluctuation in the magnitude of the drift speed 
followed by, as Wq continues to increase, a transition to 
a new state in which the direction of the drift oscillates 
periodically iFig. [231' b')'). We refer to states of this type 
as direction-reversing waves, by analogy with similar be¬ 
havior found in other systems supporting the presence 
of such waves [H, With increasing wq the rever¬ 

sals become localized in space (and possibly aperiodic, 
Fig. IMf a')') and then cease, leading to a stationary pinned 
structure at wq = 0.12 (Fig. [53Kb)). Figure [551 shows the 
position xq of the maximum of the local order parame¬ 
ter R{x, t) of the coherent state as a function of time for 
the cases in Figs. 1551 and 15^ showing the transition from 
translation to pinning as wq increases, via states that 
are reflection-symmetric on average. The final state is 
a steady reflection-symmetric pinned state aligned with 
the imposed inhomogeneity. 

In fact, the dynamics of the present system may be 
more complicated than indicated above since a small 
group of oscillators located in regions where the order 
parameter undergoes rapid variation in space may lose 
coherence in a periodic fashion even when wq = 0 thereby 
providing a competing source of periodic oscillations in 
the magnitude of the drift speed. As documented in [2^ 
this is the case when 0.7570 < /3 < 0.7644. For /3 = 0.75, 
however, the coherent state drifts uniformly when wq = 0 
and this is therefore the case studied in greatest detail. 

The profile of the pinned coherent state can also be 
determined from a self-consistency analysis. For G{x) = 
cos(x) + cos(2a::), the local order parameter Z(x) can be 
written in the form 



FIG. 23. Hidden line plot of the phase distribution ff(x, t) 
when (a) wq = 0. (b) wo = 0.04. In both cases, k = 2, 

j3 = 0.75 and N = 512. 



Z{x) = acos{x) -|-5sin(a;) -|-ccos(2a:) -b (isin(2a;). (24) 

Owing to the reflection symmetry of the solution b = d = 
0; on applying a rotation in 6 we may take a to be real. 
The self-consistency equation then becomes 

^,-0^/ cos(^)(^ + w(y)-Ml/)) \ (25) 

\ acosy + c* cos(2y) / 

_ / cos(2y)(0 -b Uj{y) - ^l{y)) \ 

\ acosy-b c* cos(2y) / 

When Wq = 0.12, k = 2 and /3 = 0.75 the solution of these 
equations is O = 1.9879, a = 3.1259, = —0.1007, Ci = 

—2.7300, results that are consistent with the values ob¬ 
tained from numerical simulation. For these parameter 
values, 11 -b w(a;) < |acos(a;) -b ccos(2a;)| for x € (—7r,7r], 
indicating that all phases rotate with the same frequency 
n. As we decrease wq to wq ~ 0.11 the profiles 11 -b oj{x) 


FIG. 24. Hidden line plot of the phase distribution 9{x, t) 
when (a) wo = 0.08. (b) wo = 0.12. In both cases, k = 2, 
P = 0.75 and N = 512. 

and |acos(x) -bccos(2a:)| start to touch (Fig. (55]) and for 
yet lower wq the stationary coherent state loses stabil¬ 
ity and begins to oscillate as described in the previous 
paragraph. 

We have also conducted experiments at a fixed value 
of Wo > 0 and k > 0 while changing /3. For example, at 
fixed Wo = 0.12, K = 2 and /3 = 0.75 the system is in 
the pinned state shown in Fig. [5^ bb Since the speed of 
the coherent state with wq = 0 gradually increases as p 
decreases, we anticipate that a given inhomogeneity will 
find it harder and harder to pin the state as /3 decreases. 
This is indeed the case, and we find that there is a critical 
value of P at which the given inhomogeneity is no longer 
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FIG. 25. The position xq of the coherent state in Figs. [23] 
and Ml as a function of time for (a) cuo = 0, (b) ojq = 0.04, 
(c) cjo = 0.08, and (d) cjo = 0.12. In all cases k = 2, /3 = 0.75 
and N = 512. 
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FIG. 26. The profiles of II+a;(a;) and i? = |a cos(a;)+ccos(2a:)| 
for ujo = 0.11, K = 2 and /I = 0.75. 



FIG. 27. Hidden line plot of the phase distribution ff(x, t) 
when (uo = 0.12 and /3 = 0.76. (a) Unpinned state obtained 
from a traveling coherent state with a;o = 0 and /? = 0.76 on 
gradually increasing oiq to 0.12. (b) Pinned state obtained 
from a traveling coherent with cuo = 0 and /? = 0.75 on grad¬ 
ually increasing ujq to 0.12, and then increasing P to 0.76. In 
both cases k = 2 and N = 512. 


(a) 


(b) 




X t 


able to pin the structure, with depinning via back and 
forth oscillations of the structure [^. For yet smaller 
values of /3, this type of oscillation also loses stability 
and evolves into near-splay states. When we increase /3 
again we uncover hysteresis in each of these transitions. 
Figure[?flshows two distinct states at identical parameter 
values: ujq — 0.12, k = 2 and /3 = 0.76 generated using 
different protocols: Fig. ETT al shows a direction-reversing 
state evolved from a traveling coherent state when we 
increase wq from 0 to 0.12 at /3 = 0.76, while Fig. HTJb) 
shows a pinned state generated from the pinned state at 
(3 = 0.75 when we change (3 from 0.75 to 0.76 at fixed 
cjq = 0 . 12 . 

B. Periodic inhomogeneity: a;(a:) = a;o cos(Za;) 

We now turn to the case uj(x) = wq cos{lx). Similar to 
the bump inhomogeneity case, when (3 is fixed, the travel¬ 
ing coherent state will continue to travel until wq reaches 
certain threshold. However, we did not find the pinned 
state as in the previous subsection around (3 = 0.75. Here 
we focus on the case (3 = 0.7 for which the traveling co- 


FIG. 28. (a) A snapshot of the phase distribution 6{x, t) of 
a coherent solution when ojo = 0.028. (b) The position a;o 
of the coherent solution as a function of time. In both cases 
1 = 1, P = 0.05 and N = 512. 


herent state has a reasonable speed when wo = 0 and 
take 1 = 1. When wq is small, the state remains coher¬ 
ent (Fig. and continues to drift, albeit no longer 

with a uniform speed of propagation. Figure[28l(b) shows 
that the speed executes slow, small amplitude oscilla¬ 
tions about a well-defined mean value v{u}q) shown in 
Fig. [29Kb); the corresponding time-averaged oscillation 
frequency H(wo) is shown in Fig. I^W ab When wq is 
increased in sufficiently small increments the oscillations 
grow in amplitude but the solution continues to travel 
to the left until ujq ~ 0.285 where a hysteretic transition 
to a near-splay state takes place. Figure ISOT bl shows 
that prior to this transition the average speed first de¬ 
creases as a consequence of the inhomogeneity, but then 
increases abruptly just before the transition owing to the 
loss of coherence on the part of a group of oscillators and 
the resulting abrupt increase in asymmetry of the order 
parameter. 
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FIG. 29. (a) The mean angular velocity and (b) the mean 
drift speed v, both as functions of a;o when I — 1, /3 = 0.7. 


V. TRAVELING CHIMERA STATES 

In addition to the states discussed in the previous 
sections, Eq. @ with constant uj admits traveling one- 
cluster chimera states when G{x) = Gn (x) = cos{nx) -h 
cos[(n -I- I)x] and appropriate values of the the phase lag 
p. This state consists of a single coherent cluster that 
drifts through an incoherent background as time evolves 
at more or less a constant speed. Figure IHOl al shows a 
snapshot of such a state when n = 3. The direction of 
motion is determined by the gradient of the phase in the 
coherent region: the cluster travels to the left when the 
gradient is positive and to the right when the gradient is 
negative. Figure ISOT bl shows the position xq of the co¬ 
herent cluster as a function of time and confirms that the 
cluster moves to the right at an almost constant speed. 
In [^ . we use numerical simulations to conclude that 
for n = 3 the traveling chimera is stable in the interval 
0.015 < /3 < 0.065. We therefore focus on the effects 
of spatial inhomogeneity on the traveling chimera state 
when P = 0.03. In fact the traveling chimera state is 
more complex than suggested in Fig. [HUf a.bl: unlike the 
states discussed in the previous sections, the profile of 
the local order parameter fluctuates in time, suggesting 
that the state does not drift strictly as a rigid object. 


A. Bump inhomogeneity: 01(2;) = ojq exp ( — K|a:|) 

In this section we investigate the effect of a bump¬ 
like inhomogeneity u}{x) = woexp (—/c|x|) on the motion 
of the traveling chimera state. Starting with the trav¬ 
eling chimera state for wq = 0, we increase uiq in steps 
of Awq = 0.01. To describe the motion of the coherent 
cluster, we follow the method in [l9l | and determine the 
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FIG. 30. (a) A snapshot of the phase distribution ff(x, t) for 
a traveling chimera state in a spatially homogeneous system, 
(b) The position xo of the coherent cluster as a function of 
time. The simulation is done for G(x) = cos(32;) -I- cos(4a;) 
with P = 0.03 and N = 512. 


instantaneous position xq of the cluster by minimizing 

N 

the function F{x*) = t) - /(xfe,x*)]^, where 

k 

f{x, X*) = — cos{x — X*) is a reference profile, and using 
the minimizer x*(t) as a proxy for xo(t). We find that 
even small ujq suffices to stop a traveling chimera from 
moving: Fig. [5T] shows that the threshold ojq ~ 0.01 for 
K = 1 and that it increases monotonically to ujq ~ 0.03 
for K = 10. The resulting pinned state persists to values 
of Wo as large as wq = 1. 


0.03 I 

0O.O2: I ^ 

o.or I J 


K 

FIG. 31. Dependence of pinning threshold oJo on k when 
P = 0.03. 

Figures |3Ua,c,e) show the position Xq of the coherent 
cluster as a function of time obtained using the above 
procedure for wq = 0.04,0.08,0.12, respectively, i.e., in 
the pinned regime. The figures show that the equilibrium 
position of the coherent region is located farther from the 
position X = 0 of the inhomogeneity peak as ujq increases. 
The bump in oj(x) thus exerts a ’’repelling force” on the 
coherent cluster, whose strength increases with the height 
of the bump. We interpret this observation as follows. 
The coherent cluster can only survive when the frequency 
gradient is sufficiently small, and is therefore repelled by 
regions where uj(x) varies rapidly. In the present case 
this implies that the coherent cluster finds it easiest to 
survive in the wings of the bump inhomogeneity, and this 
position moves further from x = 0 as wq increases. This 
interpretation is confirmed in Figs. I32f b.d.fl showing 
the average rotation frequency, 0t(x), of the oscillators. 
The plateau in the profile of 0t(x) indicates frequency 
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FIG. 32. The position xo of the pinned coherent cluster in a 
traveling chimera state as a function of time when (a) wo = 
0.04, (c) Wo = 0.08, (e) wo = 0.12. The average rotation 
frequency 6t for (b) wo = 0.04, (d) wo = 0.08, (f) wo = 0.12. 
In all cases /3 = 0.03 and k = 10, = 512. 


FIG. 33. The position xo of the pinned coherent cluster in a 
traveling chimera state as a function of time when (a) wo = 
0.005, (c) Wo = 0.01. The average rotation frequency 6t for 
(b) Wo = 0.005, (d) Wo = 0.01. In all cases /3 = 0.03, I = 1 
and N = 512. 



locking and hence the location of the coherent cluster; 
the fluctuations in the position of the coherent cluster 
are smoothed out by the time-averaging. 


FIG. 34. The position xo of the coherent cluster in a traveling 
chimera state as a function of time when (a) wo = 0.001, (b) 
Wo = 0.002. In all cases /3 = 0.03, I — 1 and N = 512. 


B. Periodic inhomogeneity: ljj { x ) = ojq cos{lx) 

We now turn to the effects of a periodic inhomogene¬ 
ity w(x) = uiocos{lx). When I = 1 and wq increases 
the coherent cluster initially travels with a non-constant 
speed but then becomes pinned in place; as in the case 
of the bump inhomogeneity quite small values of ujq suf¬ 
fice to pin the coherent cluster in place (for I = 1 the 
value Wo 0.003 suffices). As shown in Figs. [551 a. cl 
the position xq of the pinned cluster relative to the local 
maximum of the inhomogeneity (i.e., a; = 0) depends on 
the value of wq > 0.003. As shown in Figs. 1551 a. bl the 
coherent cluster travels to the right and does so with a 
speed that is larger when xo{t) > 0 than when xo{t) < 0. 
This effect becomes more pronounced as ojq increases. 
This is because the coherent structure is asymmetric, 
with a preferred direction of motion, and this asymme¬ 
try increases with wq. Evidently, the speed of the syn¬ 
chronization front at the leading edge is enhanced when 
w'(a;) < 0 but suppressed when uj'(x) > 0 and likewise 
for the desynchronization front at the rear. 

For ^ > 1 we observe similar results. The coherent 
cluster is pinned in space already at small values of wq. 


Since the inhomogeneous system has the discrete trans¬ 
lation symmetry x —>■ x + ^ the coherent cluster has I 
possible preferred positions. Figure [5S] shows an exam¬ 
ple for ^ = 2. Panels (a,c) show snapshots of the phase 
distribution 9{x,t) for wq = 0.01 in the two preferred lo¬ 
cations (separated by Aa; = tt), while panels (b,d) show 
the corresponding average rotation frequency 6t. 


VI. CONCLUSION 

In this paper we have investigated a system of non¬ 
identical phase oscillators with nonlocal coupling, focus¬ 
ing on the effects of weak spatial inhomogeneity in an 
attempt to extend earlier results on identical oscillators 
to more realistic situations. Two types of inhomogeneity 
were considered, a bump inhomogeneity in the frequency 
distribution specified by uj{x) = wq exp(—K|a;|) and a pe¬ 
riodic inhomogeneity specified by uj{x) = ujqCOs{Ix). In 
each case we examined the effect of the amplitude wq of 
the inhomogeneity and its spatial scale {l~^) on the 
properties of states known to be present in the homo¬ 
geneous case Wo = 0, viz., splay states and stationary 
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FIG. 35. (a,c) The two possible phase distributions 9{x,t) of 
pinned traveling chimera states when tuo ~ 0.01, P = 0.03 and 
I = 2. (b,d) The corresponding average rotation frequencies 
6t- In both cases N = 512. 


tude inhomogeneities sufficient to arrest the motion of 
these states. In all these cases the coherent regions are 
found in regions of least inhomogeneity, an effect that 
translates into an effective repulsive interaction between 
the coherent cluster and the inhomogeneity. 

In future work we propose to explore similar dynamics 
in systems of more realistic nonlocally coupled oscillators 
and compare the results with those for similar systems 
with a random frequency distribution. 
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Appendix A: Derivation of the self-consistency 
equation 


chimera states, traveling coherent states and traveling 
chimera states [^ . 

We have provided a fairly complete description of the 
effects of inhomogeneity on these states for the cou¬ 
pling functions G{x) = cos(^ cos(a;) -b cos(2a:) and 
cos(3a;)-|-cos(4a;) employed in [^. Specifically, we found 
that as the amplitude of the inhomogeneity increased a 
splay state turned into a near-splay state, characterized 
by a nonuniform spatial phase gradient, followed by the 
appearance of a stationary incoherent region centered 
on the location of maximum inhomogeneity amplitude. 
With further increase in ujq additional intervals of inco¬ 
herence opened up, leading to states resembling the sta¬ 
tionary multi-cluster chimera states also present in the 
homogeneous system. These transitions, like many of 
the transitions identified in this paper, could be under¬ 
stood with the help of a self-consistency analysis based 
on the Ott-Antonsen Ansatz [3l|, as described in the 
Appendix. The effect of inhomogeneity on multi-cluster 
chimera states was found to be similar: the inhomogene¬ 
ity trapped the coherent clusters in particular locations, 
and eroded their width as its amplitude ojq increased, re¬ 
sulting in coalescence of incoherent regions with increas¬ 
ing Wq. 

More significant are the effects of inhomogeneity on 
traveling coherent and traveling chimera states. Here 
the inhomogeneity predictably pins the traveling struc¬ 
tures but the details can be complex. Figures [25H23] show 
one such complex pinning transition that proceeds via an 
intermediate direction-reversing traveling wave. These 
waves are generated directly as a consequence of the in¬ 
homogeneity and would not be present otherwise, in con¬ 
trast to homogeneous systems undergoing a symmetry¬ 
breaking Hopf bifurcation as described in [s^. Many of 
the pinning transitions described here are hysteretic as 
demonstrated in Fig. [571 The traveling chimera states 
are particularly fragile in this respect, with small ampli- 


We suppose that an oscillator at position x has in¬ 
trinsic frequency oj(x), and assume that the frequency 
distribution uj(x) is continuous. The model equation is 

— = U}{x) - / G{x - y) sin(6»(x, t) - 9{y, t) + a) dy. 

(Al) 

We next introduce the probability density function 
f{x,uj,6,t) characterizing the state of the system. This 
function must satisfy the continuity equation 

df d 

^ + -(/O=0, (A2) 

where v{x,t) satisfies the relation 

/ TT 

Gix -y)f'{y,t)dy (A3) 

-TT 

and f{y, t) = sin(6> - 9' + a)f{y, w, 9', t) d9'duj. 

We also define the local order parameter 

/ TT p<X) pTT 

G{x-y) / / f{y,uj,9',t)d9'dujdy. 

-TT J —OO —TT 

(A4) 

Then 

= (A5) 

The above equations can be recast in a more convenient 
form using the Ott-Antonsen Ansatz 

/(a:,cc,0,t) = ^^[1 + + c.c.]. 

n—1 

(A6) 

Here g{x,uj) represents the distribution of natural fre¬ 
quencies at each x. Matching terms proportional to dif¬ 
ferent powers of expi9, we obtain 

^ = zcca - i [Ze-“ - , (A7) 









16 


where the complex order parameter Z(x,t) is given by 

/ TT pOO 

G{x-y) / g{x,uj)a{y,u},t)du}dy. (A8) 

-TT J —OO 


In the limit I? —> 0 the distribution fnnction g reduces 
to a delta function. The corresponding quantity z{x,t) 
satisfies 


If we take 


g{x,oj) 


D 

7r((w — uj{x)Y + D"^) 


(A9) 


dz 1 

— = iujz + -[Z exp(—iQf) — Z* exp(iQ;)2^], (All) 


set z{x,t) = a(x,uj(x) — iD,t) and perform the implied 
contour integration, we obtain 

Z{x,t)=f G{x - y)z{y,t)dy. (AlO) 

J —TT 


where Z{x,t) is given by (lAlOll . Equations (lAlOD - dAllIl 
constitute the required self-consistency description of the 
nonlocally coupled phase oscillator system with an inho¬ 
mogeneous frequency distribution uj{x). 
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